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nV«5a Project Objectives 


• Simulate Flash Thermography on Graphite/Epoxy 
Flat Bottom hole Specimen and thin void specimens. 

• Obtain Flash Thermography data on 
Graphite/Epoxy flat bottom hole specimens 

• Compare experimental results with simulation 
results 

• Compare Flat Bottom Hole Simulation with Thin 
Void Simulation to create a graph to determine size of 
IR Thermography detected defects 
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Composite Dimensions 




\ Structural Engineering Division 
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Composite Dimensions 




Structural Engineering Division^ 
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Composite Dimensions 
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Creating Flat Bottom Hole 

Simulation 



• Simulation requirements 

- Uniform thickness 

- Defects completely inside composite 

• Pixel size 

- Circular defects to square defects 

- Width in terms of pixels 


Thermal Properties 
of air 

Values 

Density 

1 .20(kg/m 3 ) 

Heat Capacity 

1 005(J/kg/K) 

Conductivity: Z axis 

0.026(W/m/K) 

Conductivity: X axis 

0.026(W/m/K) 

Conductivity: Y axis 

0.026(W/m/K) 


Thermal Properties 
of Composite 

Values 

Density 

1150(kg/m 3 ) 

Heat Capacity 

0.853(J/g/K) 

Conductivity: Z axis 

0.525(W/m/K) 

Conductivity: X axis 

3.38(W/m/K) 

Conductivity: Y axis 

3.38(W/m/K) 
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Simulation Dimensions: 

Column 1 



Specimen 

Length, L x , [m] 

0.1115 

Width, 4 [m] 

0.064 

Heat exchange coef. 
front surface, h R 
[W/(m 2 .°C)] 

10 

Heat exchange coef. 
rear surface, h R , 
[W/(m 2 .°C)] 

10 

Steps along X 

223 

Steps along Y 

128 

Number of layers, / 

1 

Number of defects 

4 

Length of each step in X, 
[m] 

5.000E-04 

Length of each step in Y, 
M 

5.000E-04 

Total thickness, L z , [in]* 

0.000 


Layers 

.ayer #1 

Conductivity, K x , [W/(m.°C)] 

0.023557 

Conductivity, K y , [W/(m.°C)] 

0.023557 

Conductivity, K z , [W/(m.°C)] 

0.0043267 

Heat capacity, C, [ J/(kg.K )] 

870.8544 

Density, r, [kg/m 3 ] 

1576.2045 

Thickness, L z m 

1 .1 18E-03 

Number of steps along Z, n 

22 

Thickness of each step in Z, [m] 

5.080E-05 

Thickness of each step in Z, [in]* 

0.002 

Timing 


Type 

Square Pulse 


Heat time, r h , [s] 

0.005 


End time, [s] 

6 


Time step, [s] 

0.005 



Heat Source 

Source in space 

Exponential 

Max heat pulse, Q, 
[W/m 2 ] 

1.800E+06 

Ambient temperature, T, 
[°C] 

30 

Initial temperature, T ( , [°C] 

30 

Coef. of spatial distribution 
in X, [1/m 2 ] 

0 

Coef. of spatial distribution 
in Y, [1/m 2 ] 

0 

Heat source center in X, [m] 

0 

Heat source center in Y, [m] 

0 

Output 

Output time step, [s] 

0.01 

Surface 

Front 
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Simulation Dimensions: 

Column 1 



Thermal Properties of Defects 

Defect A 

Defect B 

Defect C 

Defect D 

End of 
Part 

Conductivity, K x , [W/(m.°C)] 

0.026 

0.026 

0.026 

0.026 


Conductivity, K y [W/(m.°C)] 

0.026 

0.026 

0.026 

0.026 


Conductivity, K z , [W/(m.°C)] 

0.026 

0.026 

0.026 

0.026 


Heat capacity, C, [J/kg.K] 

1005 

1005 

1005 

1005 


Density, r, [kg/m 3 ] 

1.20 

1.20 

1.20 

1.20 


Length, L x , [m] 

5.500E-03 

5.500E-03 

5.500E-03 

5.500E-03 


X initial point, [m] 

2.950E-02 

5.000E-02 

7.050E-02 

9.100E-02 

1.115E-01 

Width, L y [m] 

5.500E-03 

5.500E-03 

5.500E-03 

5.500E-03 


Y initial point, [m] 

2.950E-02 

2.950E-02 

2.950E-02 

2.950E-02 


Thickness, L z , [m] 

6.096E-04 

7.620E-04 

7.620E-04 

7.620E-04 


Z initial point, [m] 

5.080E-04 

3.556E-04 

3.556E-04 

3.556E-04 
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Flat Bottom I 

Col 



• Simulation size 


D 


Front view: 



Top view: 
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Data from Flat Bottom Hole 
Simulation: Column 1 


Infrared Thermography Simulation of 0.044 Inch 
Thick Graphite/Epoxy Composite 


Max 2.751 K 


I 


/ 

Center of 
Defect D 


Center of 
Defect C 


r 

Center of 
Defect B 


Center of 
defect A 


Min 2.241 K 
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Simulation Results: 
Temperature v. Time Image 



Structural Engineering Divisions 


• Shows the difference 

in temperature. 

• Blue curve is reference 
point 



Temperature Rise (K) 


X, Y : Coordinates 
defect's center 
on simulation 



Temperature Rise 
for Reference point 


50 1 □□ 1 50 200 250 


Number of output time steps 
(Output time step 0.010000s) 
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Normalizing Data and Graphing 



Collect data from 
Temperature v. Time 
graph 

Convert text file to excel 
spreadsheet 

Normalized contrast: 

(Ti-Tj 0 )-(T r -T rn ) 

(Tj-Tj ) + (T r -T r J T = Temperature for simulation. 


Normalized Contrast: Defect A1 



O L-O LO LO t — I LO 1-0 LO f\l LO LO LO 00 LO LO LO L-O LO LO LO LO LO LO 


Time (s) 


0 


Pixel intensity for experimental IR data 
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Experimental Set-up 




Structural Engineering Division 
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Frames of Infrared 
Thermography Evaluation 




Defects A-D in Columns 1-3 
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Experimental Data 




Structural Engineering Division 


Image Window: Flat Bottom Hole 

• Reference point 

• Point of Interest 

• Different sizes 
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Simple Contrast 


Simple 



Flat Bottom Hole Contrast Evolution 



Peak Time Time (s) 
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Structural Engineering Division^ 


• Finding maximum 
simple contrast 

• Saving data as text 
file 

• Transporting data to 
Excel 

• Creating Normalized 
contrast 


7439 
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Calculating and Graphing 
Normalized Contrast 



Structural Engineering Division's/ 


• Average pre-flash 
temperatures for both 
reference point and 
point of interest 

• Use averages as 
initial temperature 


Time 

Reference 

Point 

Point of 
Interest 




-0.15 

7240 

7249 




-0.133 

7241 

7253 




-0.117 

7240 

7249 




-0.1 

7240 

7249 




-0.083 

7239 

7251 




-0.067 

7240 

7254 




-0.05 

7239 

7247 




-0.033 

7241 

7255 




-0.017 

7239 

7250 





7239.889 

7250.778 







Reference 

Point* 

Point of 
Interest* 

Normalized 

Contrast 

0 

15314 

15210 

8074.111 

7959.222 

-0.00717 

0.017 

12669 

12572 

5429.111 

5321.222 

-0.01004 

0.033 

11469 

11406 

4229.111 

4155.222 

-0.00881 

0.05 

10781 

10742 

3541.111 

3491.222 

-0.00709 

0.067 

10335 

10313 

3095.111 

3062.222 

-0.00534 

0.083 

10012 

9992 

2772.111 

2741.222 

-0.0056 

0.1 

9772 

9758 

2532.111 

2507.222 

-0.00494 
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Comparison and Correction 

Simulation 




o^r^TH^oor\imcr>c\i<£>r\iror^r^<si-oo<HLnoor\ivocr>ror s 'Cr>^i-r^ 

ooiDLfmHoootDLnmHOoocninm(NOooMnno(NOOMn 
rNjLnoq^^i^a^tNu^oq^^^^rvju^oq^m^a^rsjLnoq^rouD 
OOO^^Hr-i^-if\irsirsi roro ^^^uSuSuSuS<£><£><X) 

Time (s) 


The Normalized 
contrast for the 
original simulation 
and experimental 
data 

Analysis 

Peak Contrast 

Peak Time 

Correction of 
simulation data to 
more accurately 
portray 

experimental data 
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Correction of Simulation 




• Diffusivity: a = k /(p x C) 


• Change properties of material 


- Change Specific Heat 

- Change in Conductivity 

- Could change Density 
• Final Decision 

Simulation Data: Defect B1 

n no 

u.uo 
0.06 - 

TO °- 04 - 

C 0.02 - 

o 

U 0 - 

■O 

<u 

N -0.02 - 

Heat Capacity 

/ times 0.05 

/ \ 

— i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — — i — i — i — i — i — i — i — i 

Thermal Properties 
of Composite 

Original 

Values 

Final Values 

OfN^lDOOHfM^lDOOfNfN^ \p 00 OO (N 

OOOO^— J*— i t — I t — It— i ^ cm cm cm ^ no 

dodo dodo do c\ o o 

Density 

1150(kg/m 3 ) 

1150(kg/m 3 ) 

g -0.04 - 

^a 

o - 0.06 - 

\ 

Heat Capacity 

0.853(J/g/K) 

0.853(J/g/K) 

Z 

-0.08 - 

\ 

Conductivity: Z axis 

0.525(W/m/K) 

1.28(W/m/K) 

-0.1 - 

\ 

Conductivity: X axis 

3.38(W/m/K) 

3.85(W/m/K) 

-0.12 - 


Conductivity: Y axis 

3.38(W/m/K) 

3.85(W/m/K) 


Time (s) 
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Corrected Simulation 



Normalized Contrast: Defect A1 



OrOTH^H^HtHrsirMrMrooororo^^-^-LnLOLnLn^D<x><X)r^r^r^ 

^OOOOOOOOOOOOOOOOOOOOOOOO 

^^^(NLnoOTH^i^cDro^a^rvjLnixjTH^r^cDoo^crjrvjLn 

00^^^(N(Nrsifr)rorofr)^-^-vfLOLnLn<X)<X)<xj<X)i s ^r^ 

Time (s) 



Corrected simulation 


Comparison between 
original simulation, 
corrected simulation, 
and experimental data 
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Sources of Differences 




• Simulation contrast is based on temperature versus time. 
Experimental contrast is based on pixel intensity versus time. 

• Experimental Flash vs. Simulation Flash 

- Experimental flash envelope has a sharp rise and slow decay 

- Simulation flash is a square pulse 

• Experimental factors 

- Experimental data is more sensitive to pixel size. Get smaller pixel 
intensity for a larger pixel 

- Uneven flash causes some lateral heat flow 

- Part has a surface texture causing lateral heat flow 

• Emissivity 

- The specimen emissivity was measured to be 0.9 and provides 
lower (< 5%) experimental contrast 

• Simulation inaccuracies (model approximations, boundary condition 
approximations, no lateral heat flow) 
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Creating Thin Delaminations 



• Change depths of the defects, but leave the initial points unchanged. 

• Input data into ThermoCalc-6L 


• Run simulation 

Flat Bottom Hole Simulation: 
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Collecting Data 




Structural Engineering DMsionY 


• Same as for the flat bottom hole simulation 

- Collect data from Temperature v. Time graph for 
each defect 

- Convert the text file to excel spreadsheet 
compatible 

- Generate normalized contrast graph 
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Comparison of Simulations 



Structural Engineering Division's/ 


• Comparing flat bottom hole simulation to thin 
delamination simulation 


• Compare and graph the peak contrast ratio and 
peak time ratio 

-Thin delamination/Flat bottom hole 
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Ratio of Peak Contrast 



Peak Contrast Ratio and Peak 

Time Ratio 



Structural Engineering Division 


Peak Time Ratio vs. Peak Contrast Ratio for 
Graphite/Epoxy Composite 



Ratio of Peak Times 
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Future Work 




Structural Engineering Division's/ 


• Make controlled impacts to make thin 
delaminations 

• Evaluate delaminations with Infrared 
Thermography 

• Evaluate delaminations with Ultrasonic 
Techniques 

• Section the specimen at delaminations 

• Determine actual size of delaminations 

• Compare actual results with simulated results 

• Determine accuracy of the simulation 
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Skills Acquired 




Structural Engineering DMsiony 


• Learned Thermodynamics 

- Theory and application 

• IR temperature measurement 

• Infrared Thermography NDE 
— Simulation 

- IR Experimental data acquisition and analysis 

• Eddy Current 

• Ultrasonic Testing 

• Time management 

• Work hours 

• Technical paper 
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v%* * 

^TTieCo^ 


Experiences at JSC 


Structural Engineering DMsiony 


Volunteering 


Movie Night 



log] 








MAYJ9 W 





After Graduation 



Capestone 


Graduation 

Project 


5/2011 


Intern 
at JSC 



Co-op 

2011-2016 


Graduate School 
2011-2016 





Professor 

Work for NASA of 

Mathematics 
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Exit Presentation: Infrared 
Thermography on Graphite/Epoxy 



Thank You 
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